Evolutionary Biology

Domestic dogs (Canis familiaris) display more morphological variation than any other mammal. Cranial morphology has been extensively studied, as have the relationships with function, development, genetics, veterinary medicine, and breed welfare. Postcrania remain comparatively understudied, despite well-documented breed-specific predispositions to musculoskeletal disease. Here, we apply three-dimensional landmark-free morphometrics to quantify the shape of 743 elements from 213 dogs, including the scapula, humerus, radius, ulna, pelvic girdle, femur, tibia, and fibula. We assess integration among limb elements and investigate drivers of shape variation within and between breeds. Across most breeds, limb bone shape is strikingly similar. Dachshunds, however, exhibit distinct morphology across all elements and one to two orders of magnitude greater variation than any other breed. Despite this disparity, integration remains high between all element pairs. Remarkably, we find no significant relationship between bone shape and body mass, age, or pathology, but comparison with historic specimens reveals marked changes in dachshund long bone shape over the past [~]150 years. These extreme differences are not shared by other sampled chondrodysplastic breeds, underscoring the need to understand morphological diversity beyond simple categorisation. These findings provide a quantitative framework for linking postcranial morphology with function, disease risk, and evidence-based improvements to canine welfare.

Sarita Lake, British Columbia houses a distinctive population of threespine stickleback (Gastrosteus aculeatus L.) with a phenotype characterized by unusually large individuals relative to nearby conspecifics. We tested the hypothesis that members of this population are not isometrically larger but rather exhibit variation in allometric trajectories that reflect changes in developmental timing impacting the developmental-genetic architecture of the phenotype. We used 3D geometric morphometrics to characterize the size and shape of skulls, pectoral girdles and pelvic girdles from a sample of individuals from nearby freshwater and marine populations and compare them to a sample from Sarita Lake. We showed that individuals from the Sarita Lake population are larger in each body region compared to most other populations examined. Further, these individuals have dorsally expanded skulls and relatively robust pelvic armour. We also showed that the relationship between size and shape is differently structured among body regions and is heavily influenced by non-uniform sexually-mediated variation across populations sampled. Our results reflect complex underlying developmental trajectories, and we suggest that the large phenotype observed may be driven by fecundity selection on female size in combination with a limnetic trophic niche and relatively increased predation pressure in Sarita Lake.

How evolutionary and developmental processes interact to determine axes of neural variation that produce behavioural diversity has been debated for many decades, with alternative hypotheses giving differential emphasis to functional coupling, which favours co-evolution, and developmental constraint, which enforces it. A critical omission is data on the genetic architecture of brain size and structure, which more closely illuminates the shared developmental dependencies between components of an integrated system. Here, we exploit ecological divergence between Astatotilapia calliptera and Aulonocara stuartgranti, two closely related cichlid species from Lake Malawi, to explore the genetic architecture of brain evolution. Using computer vision and machine learning techniques to extract volumetric data from micro-tomographic images, we first demonstrate significant divergence in brain composition between these species. Genomic and micro-tomographic imaging data from a population of hybrids generated between the two species were used to investigate genetic factors shaping this differentiation. We show that the majority of brain components are integrated phenotypically in hybrids, but genetic correlations between them are generally weaker. We further show that variation in multiple brain components is associated with variation in largely structure-specific quantitative trait loci, rather than determined by genetic factors with broad effects across the entire brain. These results suggest a genetic architecture that can facilitate modular changes in brain structure, and imply that individual components are independently evolvable.

We present the first radiographic ageing framework for common dolphins (Delphinus delphis), based on ossification and epiphyseal fusion patterns in the pectoral flipper, demonstrating higher reliability for chronological age estimation than currently available epigenetic approaches for this species. Using individuals of known dental age, we calibrated two modelling approaches to predict dental age from radiographic bone scores: 1) a univariate polynomial regression using a total bone score (sum of 16 scores across all assessed flipper bones), and 2) a multivariate canonical analysis of principal coordinates (CAP) incorporating 16 individual bone-score variables. Both approaches successfully predicted dental age from skeletal ossification patterns. For an age range of 0 to 24 years, polynomial regression demonstrated high predictive accuracy with median absolute errors (MAEs) of 1.25 years in females (Spearmans {rho} = 0.93, R{superscript 2} = 0.90) and 1.08 years in males ({rho} = 0.95, R{superscript 2} = 0.86). The CAP model yielded MAEs of 1.35 years in females ({rho} = 0.90, R{superscript 2} = 0.85) and 1.80 years in males ({rho} = 0.94, R{superscript 2} = 0.84). Notably, both radiographic bone ageing models achieved equal or lower median absolute errors and higher coefficients of determination than a recently developed epigenetic clock for common dolphins derived from the same population (MAE = 1.80, Pearsons correlation (r) = 0.91, R{superscript 2} = 0.82). When applying the bone ageing models to individuals of unknown dental age, both models produced age estimates consistent with expected life-history stages (foetus, neonate, juvenile, subadult, adult), although accuracy declined in dolphins above 20 years, likely as a consequence of subtle age-related variation in skeletal changes in this species. Radiographic ageing provides an accurate non-invasive tool for demographic assessment to support conservation management of common dolphins.

Alleles with opposing effects on fitness characters are said to exhibit selectional antagonistic pleiotropy (broadly construed so that effects are not necessarily confined to the same individual). A number of theoretical investigations considered the case where a pair of alleles at a locus influences two fitness components and derived the conditions giving rise to stable polymorphism under various assumptions about the mode of trait-interaction. Strikingly, many of these analyses concluded that the potential for maintaining polymorphism is strongly constrained by the joint influence of two factors: (1) the prevalence of weak selection coefficients over coefficients of large magnitude, and (2) the absence of beneficial dominance reversals (where the deleterious effects of each allele are partially or completely masked in the heterozygous genotype). Consequently, the conclusion that selective polymorphism is unlikely to be maintained by intralocus mechanisms of antagonistic pleiotropy has achieved widespread acceptance. Here we argue that such conclusions do not apply to any of the following models of antagonism: (i) additive trait-interaction, (ii) multiplicative trait-interaction, (iii) bivoltine selection, (iv) soft selection, (v) hard selection, and (vi) sexual antagonism. We demonstrate that the parameter space giving rise to stable allelic variation is quite large throughout, and moreover, the plenitude of suitable parameters neither depends on the strength of selection nor requires dominance reversal. Dominance coefficients associated with stringent conditions for stable polymorphism are shown to be atypical as compared to all feasible parameters, and best regarded as an outcome of adherence to a special relation: dominance with a constant magnitude and direction, which includes the case of additive allelic effects at a locus. Properties of single-locus equilibria (heterozygosity, allele frequency differentiation) are investigated, as well as the contribution of dominance schemes to the genetic variance in fitness characters in populations at multilocus linkage equilibrium. Author summaryAllelic variants at a locus with opposing effects on multiple fitness components (antagonistic fitness pleiotropy) have long been appreciated as a possible source of balancing selection. The prevalence of polymorphism owing to this form of natural selection, however, has been doubted on theoretical grounds due to the fact that standard assumptions of genetic models (namely, constant magnitudes for the dominance coefficients) are hardly conducive to the maintenance of polymorphism. The major exception to this conclusion lies with schemes that exhibit dominance reversal (where the direction of dominance for antagonistic alleles flips across fitness components). Here we conduct a geometric analysis of the space of polymorphism-promoting dominance parameters and conclude that the conditions for maintaining balanced alleles is unrestrictive, with non-reversals playing an underappreciated role.

Understanding behavioral differences between non-native and closely related endangered species could be important to aid conservation management. In volume 169 of Zoology, Bels et al. (2025) reported on their comparison of display-action-patterns (DAP) between native Iguana delicatissima and non-native iguanas present on islands of the Guadeloupe Archipelago in the Caribbean Lesser Antilles. Here, we address conceptual and methodological concerns about their work and reanalyze their data given our proposed corrections, primarily a literature-informed adjustment of their "species" category. We additionally utilize online videos from South American mainland I. iguana populations, from where the non-native iguanas in the Guadeloupe Archipelago originate, to better understand the different DAPs between native and non-native iguanas in the Guadeloupe Archipelago. Significant differences in DAP characteristics among "species" categories (native I. delicatissima, non-native iguanas, and hybrids) show that Bels et al. (2025) oversimplified their data analyses by merging all non-native populations into one group. This result indicates the presence of behavioral variation among subpopulations within widely hybridizing iguanid populations, which has been poorly studied. Additionally, videos from mainland populations across two major mitochondrial clades of Iguana iguana show that non-native iguanas on Guadeloupe retained DAP characteristics of those populations from which they originate. We discuss these findings in light of the proposed hypotheses put forward by Bels et al. (2025), of which two can be excluded. Overall, our reanalysis shows that studies focusing on characteristics within settings of complex hybridization in diverse species should acknowledge this complexity.

Humans order numbers in space from left to right, with smaller quantities represented preferentially in the left hemispace and larger ones in the right hemispace. The direction of this mental number line (MNL), or more generally of number-space associations (NSA), is influenced by cultural habits such as reading and writing direction. However, a growing body of evidence from pre-verbal infants and non-human animals suggests that number-space mappings may also have biological foundations. In non-human primates, evidence for a directional MNL remains mixed, partly due to small sample sizes and methodological heterogeneity. Here, we tested samples of rhesus (Macaca mulatta) and crab-eating macaques (Macaca fascicularis) across two experiments using spontaneous food-related tasks. In Experiment 1, monkeys chose between identical food quantities (1x1 to 24x24) presented on the left and right. No systematic spatial choice bias emerged as a function of numerical magnitude, and hand use did not differ across exact numerical pairs, although exploratory analyses revealed magnitude-related modulations of manual responses. In Experiment 2, monkeys were habituated to small (4x4) or large (16x16) quantities and subsequently tested with the alternative quantity. Result showed significantly more leftward choices following numerical decreases (16[->]4) and more rightward choices following numerical increases (4[->]16), indicating that relative numerical context, rather than absolute magnitude, elicited directional spatial biases. These findings suggest that in macaques, number-space associations emerge most robustly in comparative contexts involving expectancy violations of magnitude.

At the fundamental conceptual level, two alternatives have traditionally been considered for how mutations arise and how evolution happens: 1) random mutation and natural selection, and 2) Lamarckism. Recently, the theory of Interaction-based Evolution (IBE) has been proposed, according to which mutations are neither random nor Lamarckian, but are influenced by information accumulating internally in the genome over generations. Based on the estimation-of-distribution algorithms framework, we present a simulation model that demonstrates nonrandom, non-Lamarckian mutation concretely while capturing indirectly several aspects of IBE: selection, recombination, and nonrandom, non-Lamarckian mutation interact in a complementary fashion; evolution is driven by the interaction of parsimony and fit; and random bits do not directly encode improvement but enable generalization by the manner in which they connect with the rest of the evolutionary process. Connections are drawn to Darwins observations that changed conditions increase the rate of production of heritable variation; to the causes of bell-shaped distributions of traits and how these distributions respond to selection; and to computational learning theory, where analogizing evolution to learning in accord with IBE casts individuals as examples and places the learned hypothesis at the population level. The model highlights the importance of incorporating internal integration of information through heritable change in both evolutionary theory and evolutionary computation.

Phenotypic diversity is often thought to arise from the evolutionary modification of developmental processes. However, developmental processes are tightly coupled in space and time, with each process beginning from conditions set by the one before it. While we know from dynamical systems theory that initial conditions can significantly affect a systems out-come, their importance as a source of phenotypic evolvability has been largely overlooked. Here we show for the first time, that phenotypic evolution can proceed through changes in developmental initial conditions while the underlying developmental process remains conserved. Somitogenesis is the process by which vertebral precursors, known as somites, are periodically patterned in the pre-somitic mesoderm (PSM). Somitic count (total number of somites) is thought to diversify through the evolution of components of somitogenesis such as the tempo of the segmentation clock or the mechanisms driving axial morphogenesis. Using two closely related species of Lake Malawi cichlid fishes that differ in vertebral counts, we show that somite count evolution has happened without changes to somitogenesis itself, but instead, by altering the size of the PSM at the onset of this process. This work will expand what we consider developmental drivers of phenotypic evolution and highlight the importance of comparative studies to understand the diversification of phenotypes.

We introduce a spectral existence criterion for the evolution of cooperation in the form of the inequality{lambda} maxb > c, where{lambda} max is the leading eigenvalue of an interaction operator encoding population structure, and b and c represent benefit and cost tradeoffs, respectively. Nowaks five rules for the evolution of cooperation correspond to cases in which the cooperation condition reduces to a scalar assortment coefficient. These results follow from the Price equation, which sheds light on a long-standing debate on the role of inclusive fitness and evolutionary dynamics in explaining the evolution of cooperation.

Understanding the distribution of paternity within social groups is critical for testing hypotheses about the evolution of behavior and morphology in primates, but assembling the requisite comparative data is a challenging task. We compiled genetic paternity data from 52 species of wild nonhuman primates along with information about socioecological, morphological, and life history traits that are relevant to understanding what proportion of offspring are sired by primary males (i.e., alpha males in multi-male groups and resident males in single male groups). Our dataset, which currently contains information about 11 primate families and >3,000 individual paternities, is presented as a publicly accessible, living database designed to be updated as new data become available. Using Bayesian regression models, we investigated the role that phylogeny, group composition, and seasonality play in determining primary males paternity share, and assessed the relative share of paternities obtained by non-primary residents versus extra-group males. First, we found that phylogeny has a detectable but relatively modest influence on primary males paternity share. Species-level differences explained roughly 35-40% of variation in primary males paternity share, and of that interspecific variation, [~]50-70% was attributable to shared phylogenetic history. Second, group composition strongly predicted paternity share outcomes. Primary males in single-male/multi-female groups obtained the highest share of paternity ([~]80%), while those in multi-male groups had the lowest ([~]60%), though there was substantial variation within each category. Pair-living animals showed a striking split: males in cohesive pairs sired [~]90% of offspring, while those in dispersed pairs sired only [~]55%. Contrary to expectations, reproductive seasonality did not predict primary males paternity share in any group type. Finally, when primary males in multi-male groups lost paternities, [~]75% of losses were to other resident males. Overall, [~]5-15% of offspring in these groups were sired by extra-group males. Our results largely confirm earlier findings based on smaller datasets, but also show that the relationship between social organization and paternity is more complicated than simple categorical predictions suggest. We discuss the gap between the data that would ideally be available for testing these hypotheses versus what currently exists, with hopes that our living database can help close this gap over time.

Our bodies have evolved to maintain homeostasis through regulatory systems that continuously adapt to keep physiological processes within a normal range. From this perspective, complex chronic disease can be understood as a breakdown of compensatory mechanisms, resulting in loss of homeostasis. Here we propose that adaptive receptor expression dynamics may serve as one such compensatory mechanism, increasing receptor surface expression when external ligand is insufficient, and clearing it when signaling is excessive. To explore this, we adapt a previously published agent-based model and simulate it under a range of scenarios. We find that the system of adaptive receptor expression is robust to oscillatory perturbations but not to chronic stress. We propose that receptor turnover dynamics may be better understood as an adaptive, environmentally responsive process rather than a fixed biological property, and that in some cases, disease manifests only after compensatory mechanisms have been pushed past their limits. We conclude with a discussion of implications for understanding complex chronic diseases, for thinking about epigenetic and mutational change as escalating layers of adaptation, and for how we model receptor dynamics in the context of receptor-mediated drug activity.

Parental rejection of apparently healthy newborns is widely classified as a behavioural abnormality in captive primate colonies, yet its biological significance remains unclear. In owl monkeys (Aotus nancymaae), parental rejection, defined here as cessation of nursing leading to rescue nursery rearing, is typically lethal for offspring and is transmitted across generations despite reducing offspring survival. Here, we tested whether parental rejection is associated with lifespan and reproductive differences in parents and their surviving offspring. We analysed long-term demographic records from a captive colony of 962 individuals and compared survival and reproductive outcomes between rejector and non-rejector parents using survival analyses and regression-based models. Parents who rejected offspring lived significantly longer than non-rejectors, with an average lifespan advantage of approximately 4-4.5 years in both males and females. This survival difference was concentrated during the prime reproductive period (6-20 years). Well-reared offspring of rejector parents also lived longer than offspring of non-rejectors, with a mean lifespan difference of 1.26 years. Rejector parents produced more offspring overall, but this difference was explained by extended lifespan rather than higher reproductive output per year. Analyses stratified by rejection timing showed no longevity advantage in first-birth rejectors, whereas parents rejecting later-born offspring exhibited longer survival. Together, these findings show that parental rejection is associated with longer lifespan in parents and in their well-reared offspring under captive conditions. These patterns are consistent with altered allocation of parental investment under energetic or environmental stress.

Numerical abilities are widespread in the animal kingdom and are not exclusive to humans. Domestic chicks (Gallus gallus) have been shown to discriminate numerosities spontaneously, but prior research has focused exclusively on the visual modality. Whether chicks can discriminate numerical information in the auditory domain remains unknown, despite evidence that they can perceive other auditory features such as tone and rhythm. In this study, we investigated spontaneous numerical discrimination in the auditory modality in naive domestic chicks. In Experiment 1, newly-hatched chicks were tested for their ability to discriminate between two auditory sequences differing in numerosity (4 vs. 12 identical sounds), with and without controlling for continuous variables such as duration and total sound amount. Experiment 2 examined chicks filial imprinting responses to familiar or unfamiliar numerosities. Experiment 3 controlled for potential spontaneous preferences for a single longer sound versus a shorter one. Our results showed a preference for the 12-sound sequence only when duration and total sound amount were not matched. When these continuous variables were controlled, no spontaneous numerical preference emerged. Experiment 2 revealed an overall preference for the 12-sound sequence regardless of imprinting conditions, while Experiment 3 confirmed that chicks do not have an inherent preference for longer sounds. These findings suggest that chicks are sensitive to overall magnitude in the auditory domain but do not spontaneously discriminate numerical differences when other continuous variables are held constant. Future studies will explore how specific stimulus features, such as heterogeneity of sounds, influence these preferences.

According to sexual selection theory, males should benefit more from mating with multiple partners than females do, as male investment into offspring production is typically lower. For females, empirical evidence indeed often shows diminishing returns or even costs of mating multiply. For males, the assumption often seems to be "the more, the better" - i.e., a steady increase of male reproductive success with mate number - but experimental tests of it are rare. Here we used a laboratory experiment with Trinidadian guppies (Poecilia reticulata), known for being promiscuous, to assess how pairing males weekly with 4 vs. 7 females affects both sexes reproductive performance (n = 32 polygynous males and 170 monogamous females). Increased polygyny delayed females reproductive onset by 9% and tripled their risk of reproductive failure. High-polygyny males fathered offspring with 49% more females and had 73% higher daily reproductive output. Yet, they needed 19% longer to initiate pregnancy, and only accumulated more offspring than low-polygyny males after two months. This study suggests that male mating performance is not unlimited. Especially when high extrinsic mortality selects for fast reproduction, less polygyny might be advantageous, and the strength of sexual selection perhaps more similar between the sexes than often assumed.

Dermal bone, which forms a variety of skeletal structures and persists in a wide range of extant vertebrates, evolved prior to endochondral bone which forms all mammalian load-bearing bones. Sturgeons are a family of fish which diverged soon after the lobe-finned/ray-finned split. Sturgeon retain a long robust spine at the leading edge of the pectoral fin, called the pectoral fin spine (PFS). Pectoral fin spines are bone elements that are present in many extinct and extant species of non-tetrapod jawed fish. In this study, we characterize the structure (light, polarized, micro-computed tomography and scanning electron microscopy), composition (FTIR, TGA, BMD), and mechanical properties (3-point bending and microindentation) of the pectoral fin spine (PFS) of the Russian sturgeon (Huso gueldenstaedtii). The microstructure of the PFS is highly organized as it is formed by dermal osteonal bone and parallel fibered bone. Its microarchitecture, along with high material toughness, anisotropy, and substantial ash content, enables the PFS to bear loads and function in both locomotion and protection. In addition, we show an interconnected network of neurovascular canals and ornamentations, features also found in pectoral fin spines of other non-tetrapod jawed fish. Collectively, these findings demonstrate that dermal bone can form structurally organized, mechanically competent load-bearing elements and provide new insight into pectoral fin spines in ray-finned fish.

Socially-controlled sex changing fishes provide powerful model systems for investigating sexual development and phenotypic plasticity in both behavior and physiology. The remarkable sexual transformation these fishes undertake is strongly influenced by their position in dominance hierarchies. However, the behavioral mechanisms underlying hierarchical formation remain understudied, particularly among female groups. Here, we investigated the role of winner-loser effects among females in establishing social dominance in a female-to-male sex changing fish. Individuals with prior losing experiences were more likely to lose subsequent size-matched fights, demonstrating clear loser effects, while there was no evidence for winner effects. Initial mirror aggression and some prior fighting behaviors, particularly submission, significantly and positively correlated with aggression in size-matched fights and subsequent mirror aggression; however, contest outcomes were not altered by these factors. Additionally, mirror aggression increased significantly only in subjects that drew size-matched fights. These findings demonstrate complex fighting dynamics in female-female competition and confirm the presence of loser effects in a sequential hermaphroditic species. These effects may represent evolutionarily advantageous mechanisms underlying sex change, thereby offering further context for examining how social rank advantages drive sexual transition.

Dental enamel, the hardest mineralised tissue in the human body, has proven to be an excellent source of ancient proteins, which have been found to survive within dental enamel for at least twenty million years. In archaeological and palaeontological contexts, the enamel proteome is generally considered to be rather small, consisting of about twelve proteins, most of which are unique to enamel. During amelogenesis these proteins undergo in vivo digestion by matrix metalloproteinase 20 (MMP20) and kallikrein 4 (KLK4) as well as serine phosphorylation by family with sequence similarity member 20-C (FAM20C) that alter their characteristics. Gaining knowledge of the previously understudied influence of amelogenesis on the archaeological human dental enamel proteome could benefit various palaeoproteomic analysis, especially in an human evolutionary context. Here we present archaeological dental enamel proteomes and explore protein cleavage patterns and sequence coverage to estimate the effects of in vivo digestion, as well as explore phosphorylation patterns. Additionally, we present a new marker based on phosphorylation to estimate genetic sex.

Evolving populations, especially in the strong-selection-weak-mutation limit, can be modelled as adaptive walks on fitness landscapes, moving in fitness-increasing mutational steps until reaching a fitness peak--a local optimum. Simulations of such adaptive walks--on a multi-peaked empirical landscape of the folA gene and on landscapes generated by the Rough Mount Fuji (RMF) model-- have shown that some landscapes are highly navigable, meaning that the highest x% of peaks are reached by >> x% of adaptive walks. This prompts the question of how adaptive walks can be so successful despite the local, myopic rules behind each adaptive step. Here, we investigate this question using simulations and mathematical approximations of random adaptive walks on a simplified RMF landscape. The landscape has a low-to-intermediate fitness region, whose size reconciles a low peak density with a high peak number. Despite the high number of peaks, walkers are likely to exit this region without terminating at a peak because the probability of a peak transition at each step is low and a fitness gradient guides walkers to the high-fitness region in few steps. Thus, three features are sufficient to explain why adaptive walks in the simplified RMF landscape are likely to reach a small fraction of top-ranking peaks: a low-to-intermediate fitness region with a high number of peaks, a low peak-transition probability, and which is crossed in few steps. We find that these three features are also present in the empirical folA landscape, suggesting that similar principles may apply.

Life history traits are often correlated, creating trade-offs that may impede the response to natural selection and be responsible for the evolution of senescence. These trade-offs may arise through pleiotropic effects, which can affect the response to selection in ways that resemble intra-locus sexual antagonism. Despite these hypothesized relationships, we lack clear connections between pleiotropy, sexual antagonism, and the evolution of life histories. Empirical tests for inter-sexual differences in life-history traits, including sex-specific aging, can be used to evaluate hypotheses about how pleiotropy and sexual conflict affect evolutionary trade-offs. To those ends, we measured lifespan, development time, and body size in Drosophila pseudoobscura males and females, each of which carried one of six third chromosome inversion genotypes. Temperature affected lifespan and development more than any other factor; higher temperatures increased mortality rate, decreased lifespan, and accelerated development. However, we also observed sex differences in mortality rates and development times that depended on genotype and temperature. Notably, temperature elevated the initial mortality rate across all flies, yet increasing temperatures reduced the rate of aging in some genotype-sex combinations. Similarly, direct effects of genotype on mortality rate and development time depended greatly on sex and temperature, but there was no genotype effect on body size. Despite these context-dependent genotype effects on life history traits, we failed to identify any correlations that would serve as clear evidence for sexual conflict or trade-offs. Our results therefore suggest that either historical conflicts have been resolved or any conflicts that may exist do not result in the correlations predicted by existing models.